Joint constraints on the lepton asymmetry of the Universe and neutrino mass from 

the Wilkinson Microwave Anisotropy Probe 
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We use the Wilkinson Microwave Anisotropy Probe (WMAP) data on the spectrum of cosmic 
microwave background anisotropics to put constraints on the present amount of lepton asymmetry 
L, parameterized by the dimensionless chemical potential (also called degeneracy parameter) £ and 
on the effective number of relativistic particle species. We assume a flat cosmological model with 
three thermally distributed neutrino species having all the same mass and chemical potential, plus an 
additional amount of effectively massless exotic particle species. The extra energy density associated 
to these species is parameterized through an effective number of additional species AN°^ CIB . We 
find that < |£| < 1.1 and correspondingly < \L\ < 0.9 at 2a, so that WMAP data alone cannot 
firmly rule out scenarios with a large lepton number; moreover, a small preference for this kind 
of scenarios is actually found. We also discuss the effect of the asymmetry on the estimation of 
other parameters and in particular of the neutrino mass. In the case of perfect lepton symmetry, we 
obtain the standard results. When the amount of asymmetry is left free, we find m v < 3.6 eV at 
2a. Finally we study how the determination of \L\ is affected by the assumptions on AN^ eTS . We 
find that lower values of the extra energy density allow for larger values of the lepton asymmetry, 
effectively ruling out, at 2a level, lepton symmetric models with AA r °^ lors ~ 0. 

PACS numbers: 98.70.Vc, 14.60.Pq, 98.80.Es 



I. INTRODUCTION 



It is a remarkable fact that our observational knowl- 
edge of the Universe can be justified in terms of a model, 
the so-called power-law ACDM model, characterized by 
just six parameters, describing the matter content of our 
Universe (the physical density of baryons u>b, the physical 
density of matter u> m , the Hubble constant h), the initial 
conditions from which it evolved (the amplitude A and 
the spectral index n of the primordial power spectrum) 
and the optical depth at reionization (t). In particu- 
lar, this model provides a good fit to both the cosmic 
microwave background (CMB) [l| and large scale struc- 
ture (LSS) data (although in this last case one additional 
parameter, the bias parameter b, is needed) 0. Nev- 
ertheless, the data leave room for more refined models, 
described by additional parameters: among them, the 
spatial curvature, the amplitude of tensor fluctuations, a 
running spectral index for scalar modes, the equation of 
state for dark energy, the neutrino fraction in the dark 
matter component, a non-standard value for the relativis- 
tic energy density. All have been considered in previous 
works. In particular the last two have been studied in or- 
der to gain deeper information on the properties of neu- 

trmos QiifiiaaiiamEiaGiiaeiia. 

Before the measurements of the CMB anisotropy spec- 
trum carried out by the Wilkinson Microwave Anisotropy 
Probe (WMAP) [J E E HI WL El E|, the com- 
bined CMB and LSS data yielded the following upper 



bound on the sum of neutrino masses: m v ^ 3 eV 
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The WMAP precision data allowed to strengthen this 
limit. Using rather simplifying assumptions, i.e., assum- 
ing three thermalized neutrino families all with the same 
mass and a null chemical potential (thus implying perfect 
lepton symmetry), the WMAP team found that the neu- 
trino mass should be lower than 0.23 eV |l|. This tight 
limit has been somewhat relaxed to Y m u < 1 eV Q ow- 
ing to a more careful treatment of the Ly-cv data, and its 
dependence on the priors has been examined [5j . The LSS 
data can also be used to put similar constraints, although 
they are usually weaker. Using the data from the 2 De- 
gree Field Galaxy Redshift Survey (2dFGRS) and assum- 
ing "concordance" values for the matter density Q m and 
the Hubble constant h it is found that J2m u < 1.8 eV 
0. A combined analysis of the Sloan Digital Sky Sur- 
vey (SDSS) and WMAP data gives a similar bound : 
^m„ < 1.7 eV Q ■ Quite interestingly, the authors of 
Ref. [2j claim that, from a conservative point of view (i.e., 
making as few assumptions as possible) , the WMAP data 
alone don't give any information about the neutrino mass 
and are indeed consistent with neutrinos making up the 
100% of dark matter. In Ref. it is claimed that the 
cosmological data favor a non-zero neutrino mass at the 
68% confidence limit, while the authors of Ref. [8( find 
the limit J2 m v < °- 74 - 

At the same time, more detailed scenarios with a dif- 
ferent structure of the neutrino sector have been stud- 
ied. The first and more natural extension to the stan- 
dard scenario is the one in which a certain degree of 
lepton asymmetry (parameterized by the so-called de- 
generacy parameter £, i.e., the dimensionless chemical 
potential) is introduced [2J, [2{j |26j- Although stan- 
dard models of baryogenesis (for example those based 
on SU(5) grand unification models) predict the lepton 
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charge asymmetry to be of the same order of the baryonic 
asymmetry B ~ 1CP 10 , nevertheless there are many par- 
ticle physics motivated scenario in which a lepton asym- 
metry much larger than the baryionic one is generated 
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cases, the predicted lepton asymmetry can be of order 
unity. One of the interesting cosmological implications 
of a net leptonic asymmetry is the possibility to gener- 
ate small observed baryonic asymmetry of the Universe 
|40ll4l| via the so-called sphaleron process 0]- The pro- 
cess of Big Bang Nucleosynthesis (BBN) is very sensitive 
to a lepton asymmetry in the electronic sector, since an 
excess (deficit) of electron neutrinos with respect to their 
antiparticles, alters the equilibrium of beta reactions and 
leads to a lower (higher) cosmological neutron to proton 
ratio n/p. On the other hand, an asymmetry in the /j, 
or r sector, even if not influencing directly the beta re- 
actions, can increase the equilibrium n/p ratio due to a 
faster cosmological expansion. This can be used to con- 
strain the value of the degeneracy parameter . This 
leads to the bounds -0.01 < £ e < 0.22 and |^ T | < 2.6 
011. 

The effect of a relic neutrino asymmetry on the CMB 
anisotropy and matter power spectrum was first studied 
in Ref. [46] ■ and is mainly related to the fact that a lep- 
ton asymmetry implies an energy density in relativistic 
particles larger than in the standard case. The cosmolog- 
ical observables can then be used to constrain this extra 
energy density, parameterized by the effective number of 
relativistic neutrino species N . Although this is some- 
what more general than the case of a lepton asymmetry, 
in the sense that the extra energy density can arise due to 
other effects as well, nevertheless the case of a non-null 
chemical potential is not strictly covered by the intro- 
duction of N eS . This is because the increased relativistic 
energy density is not the only effect connected to the lep- 
ton asymmetry (an additional side effect is for example 
a change in the Jeans mass of neutrinos 
the hypothesis of a negligible neutrino mass, it has been 
shown that the WMAP data constrain N cS to be smaller 
than 9; when other CMB and LSS data are taken into 
account, the bound shrinks to 1.4 < N eS < 6.8 HE3- A 
combined analysis of CMB and BBN data leads to even 
tighter bounds 0, 0, 0, 0. A more detailed anal- 
ysis, in which the effective number of relativistic relics 
and the neutrino mass are both left arbitrary and varied 
independently, can be found in Ref. ^4|. I n the same 
paper, the effect of different mass splittings is also stud- 
ied. Finally an extension of these arguments to the case 
in which additional relativistic, low-mass relics (such as 
a fourth, sterile neutrino or a QCD axion) are present, 
has been studied in Ref. |l5j . 

The goal of this paper is to perform an analysis of the 
WMAP data using the degeneracy parameter, together 
with the effective number of relativistic particles, as addi- 
tional free parameters, in order to put constraints on the 
lepton number of the Universe. We work in the frame- 
work of an extended cosmological model with three ther- 



mally distributed neutrino families having all the same 
mass and chemical potential, plus a certain amount of 
exotic particles species, considered to be effectively mass- 
less. We use the physical neutrino density lo v = £l„/i 2 , 
the degeneracy parameter £ and the extra energy den- 
sity in exotic particles AiV^ as additional parameters 
that describe the neutrino sector. We perform an analy- 
sis in a 8-dimensional parameter space that includes the 
standard, "core" cosmological parameters. 

The paper is organized as follows. After a discussion on 
the motivations that drive our work in Sec. [HI we shortly 
review some basic formulae in Sec. IIIII and discuss the 
impact of a non-null degeneracy parameter on the CMB 
spectrum in Sec. IIVI In Sec. |V|we describe the analysis 
pipeline, while in Sec. IVII we present our basic results. 
Finally, we draw our conclusions in Sec. IVIII 



II. MOTIVATION FOR THIS WORK 

The main motivation for this work comes from the fact 
that, even though several analyses have been performed 
which were aimed at putting constraints on the number 
of effective relativistic degrees of freedom, a statistical 
analysis of the CMB data aiming to put bounds directly 
on the degeneracy parameter, instead of on N , is nev- 
ertheless still missing. There are two reasons for this: 
first of all, in the limit of a vanishing neutrino mass, the 
increase in N is in effect all that is needed to imple- 
ment the non null chemical potential into the standard 
model of the evolution of perturbations [ZtiL |48| . It is 
then argued that, since neutrinos with mass smaller than 
roughly 0.3 eV, being still relativistic at the time of last 
scattering, would behave as massless, the distinction be- 
tween £ and 7V eff is no more relevant in this case for what 
concerns their effect on the CMB anisotropy spectrum. 
Although this is certainly true, it is our opinion that this 
does not allow to neglect a priori the difference between 
the two parameters. One reason is that the most conser- 
vative bound on neutrino mass, coming from the tritium 
beta decay experiments, reads m u < 2.2 eV (at the 2a 
level) |43,|53, that is quite higher than the value of 0.3 eV 
quoted above. The main evidences for a neutrino mass 
in the sub-eV range come, in the field of particle physics, 
from the experiments on neutrinoless double beta decay 
|5lL l52| . whose interpretation depends on assumptions 
about the Majorana nature of neutrinos and on the de- 
tails of the mixing matrix. Other indications of a sub-eV 
mass come, as stated above, from cosmology and in par- 
ticular from the power spectrum of anisotropies, but since 
we want to keep our results as much as possible indepen- 
dent from other analyses, we should not use information 
on neutrino mass derived from the previous analyses of 
the CMB data. Moreover, let us note that CMB data 
analyses are often refined using the results from LSS ex- 
periments. Since the structure formation, starting close 
to the epoch of matter-radiation equality, goes on until 
very late times, even very light neutrinos (in the range 
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10~ 3 -^0.3 eV) cannot be considered massless for the pur- 
pose of evaluating their effect on the matter power spec- 
trum. This means in particular that using N cS would 
lead to overlook the change in the free streaming length 
and in the Jeans mass of neutrinos due to the increased 
velocity dispersion 53] . It is then our opinion that the 
use of N eS , even if correct with respect to the interpreta- 
tion of CMB data, precludes the possibility of correctly 
implementing the LSS data as a subsequent step in the 
analysis pipeline. 

The second point against the cosmological significance 
of the degeneracy parameter is related to the constrain 
from BBN. It was recently shown that, if the Large Mix- 
ing Angle (LMA) solution to the solar neutrino problem 
is correct fas the results of the KamLAND experiment 
suggest then the flavor neutrino oscillations equal- 

ize the chemical potentials of e, /i and r neutrinos prior 
to the onset of BBN, so that a stringent limit £ < 0.07 
actually applies to all flavours [55|, |56|, I&3 • This would 
constrain the lepton asymmetry of the Universe to such 
small values that it could be safely ignored in cosmo- 
logical analyses. However, the presence of another rel- 
ativistic particle or scalar field would make these limits 
relax [5^ , while the effect of the mixing with a light ster- 
ile neutrino, whose existence is required in order to ac- 
count for the results of the Liquid Scintillation Neutrino 
Detector (LSND) experiment |5^], is still unclear |5t| . 
Moreover, it has been recently shown that a hypothetical 
neutrino-majoron coupling can suppress neutrino flavor 
oscillations, thus reopening a window for a large lepton 
asymmetry jfKjEOJ. F° r a U these reasons, we judge it is 
interesting to study if CMB data alone can constraint or 
maybe even rule out such exotic scenarios. 



III. BASIC FORMULAE 

It is customary in cosmology to call ultrarelativistic (or 
simply relativistic) a species x that decouples from the 
photon bath at a temperature Td such that its thermal 
energy is much larger than its rest mass energy: ksTd 
m x c 2 . 

Owing to Liouville's theorem, the distribution func- 
tion in momentum space f x (p ; T x , £ x ) of the species x is 
given, after decoupling, by (we shall use all throughout 
the paper units in which c = h = ks = 1): 
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where £ = /id/Td is the dimensionless chemical potential, 
often called degeneracy parameter, the sign + (— ) corre- 
sponds to the case in which the x's are fermions (bosons), 
g is the number of quantum degrees of freedom, and the 
temperature T evolves in time as the inverse of the cos- 
mological scale factor a, so that T(t) ■ a(t) = const. 
The energy density of the x's at a given temperature 



is readily calculated: 



p x (T x ,0= / E(p)f(p;T x ,Z x )d 3 p 
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p 2 ^/p 2 +mlf(p;T x ^ x )dp. (2) 



Using the dimensionless quantities y = p/T and 
f3 = m x /T, the expression for the energy density can 
be put in the form: 



9 f°° 
Px(T x ,£ x ) = -^ T x J dyy 1 



exp(y-£ x ) ± 1 



(3) 



We stress the fact that a temperature dependence is still 
present in the integral through the term [3. However, the 
temperature dependence disappear from the integral in 
two notable limits, the ultrarelativistic (UR) and non- 
relativistic (NR) one, corresponding respectively to the 
two opposite cases (3 <C 1 and [3 3> 1 J2{| . Then, defining 



±fc\ — 



y" 



e»~£ ± 1 



dy 



e v ± 1 



dy) , (4) 



so that J„(Q) = 1, we have 



Px{T x ,£ x ) = 



7/8 ) 9*h J f(^ T * UR 



(3/4) 9*^rrn x j£(t x )T* NR 



(5) 

where the upper and lower values in parentheses in front 
of the expression in the right-hand side hold for bosons 
and fermions respectively, and Ci n ) is the Riemann Zeta 
function of order n. 

It is useful to express p x (t) in terms of the present day 
energy density of the cosmic background photons: 



Px(t) = 



1 

7/8 
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having defined an effective number of relativistic degrees 
of freedom g° s as 
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It is often the case that one has to consider a fermion 
species x together with its antiparticle x, the most no- 
table example being the relic neutrinos and antineutri- 
nos. In chemical equilibrium, the relation £ x — —£ x holds 
owing to the conservation of chemical potential, as can 
be seen considering the reaction: 



x + x 



7 + 7 



(8) 
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and noting that the chemical potential in the final state 
vanishes |62| . This relation holds for neutrinos and an- 
tineutrinos in several cosmological scenarios. There are 
some exceptions to this, most notably early Universe sce- 
narios in which lepton asymmetry is generated |63l | or de- 
stroyed by active-sterile neutrino oscillations at low 
temperatures. However, we shall assume all throughout 
the paper that the relation £„ = — £ p holds. 
It can then be shown that 



followed the electron-positron annihilation, shortly after 
neutrino decoupling. Then 



9 =g (77 J N„* 0.23 N„, (13) 

where N v = 3 is the number of neutrino families. The 
energy density in a single neutrino species is: 
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where the factor between square parentheses can be rec- 
ognized as what it is often quoted as the contribution of a 
non-vanishing chemical potential to the effective number 
of relativistic species iV cff . 

The definitions introduced above can be easily ex- 
tended to the case when several ultra-relativistic species 
Xi are present: 
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where photons are excluded from the summation. This 
means that, since g° s — g 1 = 1, the actual number of 
relativistic degrees of freedom is (1 + g ). 

The total density of ultrarelativistic particles at a given 
time is thus: 



Aad = P° (1 + 5 cff ) (l + z) 4 



(11) 



Finally we can use this expression to find the depen- 
dence on g cS of the redshift of radiation-matter equality 
z eq (the subscripts b and CDM stands for baryons and 
cold dark matter respectively): 
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(14) 



However, several mechanisms that could increase (or 
even decrease) the energy density of relativistic parti- 
cles have been proposed. In the presence of some extra 
relics (such as sterile neutrinos, majorons, axions, etc.) 
the energy density of radiation would obviously increase. 
A non-zero chemical potential for neutrinos or an unac- 
counted change of p 1 , due for example to particle decays 
that increase the photon temperature, would produce the 
same result. In all cases the effect is the same: a change 
in g cS , as it can be seen by looking at eq. 0. It is usual 
in the literature to parameterize the extra energy density 
by introducing an effective number of neutrino families 
7V eff , defined as: 



N 



ctr 



n std 
rv 



(15) 



where again the sum runs over all ultrarelativistic 
species with the exceptions of photons. It is clear from 
this definition that N eS is actually the energy density in 
ultrarelativistic species (apart from photons) normalized 
to the energy density of a single neutrino species with 
zero chemical potential and standard thermal history. It 
is easy to show that a relation formally similar to (|13|) 
holds in the non-standard scenario: 



l + z eq = Pb+P Q c ™ (l + fl**)- 1 . (12) 
Pi 

So, the larger is the energy density of ultra-relativistic 
particles in the Universe, parameterized by the effective 
number of degrees of freedom g , the smaller z eq will be, 
i.e., the later the equality between radiation and matter 
will occur. In other words, supposing that the density in 
non-relativistic particles (baryons + CDM) is well known 
and fixed, having more relativistic degrees of freedom will 
shift z eq closer to us and to the CMB decoupling. 

In the standard cosmological scenario, the only contri- 
bution to the energy density of relativistic particles other 
than photons is the one due to the three families of stan- 
dard neutrinos, with zero chemical potential. The ratio 
of the neutrino temperature to the photon temperature 
is T°/T° = (4/11) 1 / 3 , due to the entropy transfer that 



g ca = 0.23iV ctt . (16) 

In addition, it should be noted that even in the stan- 
dard scenario N cS ^ N v = 3 ) but instead N cS ~ 3.04. 
This is due to the fact that neutrino decoupling is not in- 
stantaneous, so that neutrino actually share some of the 
entropy transfer of the e + e _ annihilation, on one side, 
and to finite temperature quantum electrodynamics cor- 
rections on the other |65|, |6(| . 

It is also useful to introduce the effective number of 
additional relativistic species AN cS defined as: 

AiV cff = iV cff - 3.04 , (17) 

so that AN = in the standard scenario. Please note 
that AN eS can also be negative, for example in very low 
reheating scenarios |67| . 
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In this paper we shall consider a scenario in which the 
radiation content of the Universe at the time of radiation- 
matter equality is shared among photons, three neutrino 
families with standard temperature but possibly non- 
zero chemical potential, and some other relic particle. 
We shall suppose that the presence of the latter can be 
completely taken into account through its effect on N cS . 
This is true if the species has been in its ultrarelativistic 
regime for the most part of the history of the Universe. 
The presence of this extra relic is required for our anal- 
ysis, in order to circumvent the equalization of neutrino 
chemical potentials, as explained at the end of sectioulTTl 
We also assume that the degeneracy parameters for neu- 
trinos and antineutrinos are equal and opposite, and that 
e, /x and r neutrinos all have the same chemical potential. 

The extra energy density can thus be split into two 
distinct contributions, the first due to the non-zero de- 
generacy parameter of neutrinos and the second due to 
the extra relic (s): 



AiV off = AN f(£) + AN; 



others' 



(18) 



Following our assumptions, AN° S can be expressed as a 
function of the chemical potential only: 



ANf (0 = 3 



(19) 



IV. EFFECT OF A NON-NULL CHEMICAL 
POTENTIAL 

As anticipated above, the main effect connected to the 
presence of a non- vanishing degeneracy parameter, is an 
increase in g eS (or, equivalently, in N eS ). The presence of 
this extra number of effective relativistic degrees of free- 
dom can in principle be detected from observations of 
the CMB radiation. The shift of matter-radiation equal- 
ity has important consequences for the CMB anisotropy 
spectrum, these being due to the larger amplitude of the 
oscillations that enter the horizon during the radiation 
dominated phase, and to a larger early integrated Sachs- 
Wolfe (ISW) effect. However these effects, basically due 
to the speeding up of the cosmological expansion, can be 
similarly produced by the variation of other cosmological 
parameters, for example by a smaller CDM density. 

Moreover since the change in the redshift of matter- 
radiation equality depends on the ultra-relativistic 
species only through the quantity g cS , it cannot be used 
to distinguish between the different species (i.e., it is "fla- 
vor blind"), nor to understand if the excess energy den- 
sity is due to the presence of some unconventional relic, 
to an extra entropy transfer to photons, or to a non- 
null chemical potential (i.e., to a lepton asymmetry), or 
maybe to all of the previous. 

However ultrarelativistic particles, other than chang- 
ing the background evolution, have an effect even on the 
evolution of perturbations, as it was pointed out in |6S| 



with particular regard to the case of neutrinos. First of 
all, the high velocity dispersion of ultrarelativistic par- 
ticles damps all perturbations under the horizon scale. 
Second, the anisotropic part of the neutrino stress-energy 
tensor couples with the tensor part of the metric per- 
turbations. It was shown in [£9j that this reduces the 
amplitude squared of tensor modes by roughly 30% at 
small scales. Finally, the authors of Ref. [68( claimed 
that the perturbations of relativistic neutrinos produce a 
distinctive phase shift of the CMB acoustic oscillations. 
These effects can thus be used to break the degeneracy 
between g eS and other parameters. It remains to estab- 
lish whether they can be used to break the degeneracy 
between the different contributions to g or not. Even 
without performing a detailed analysis, it can be see n by 
looking at the relevant equations in Ref. |68l | and |69| . 
that both the absorption of tensor modes and the phase 
shift depend on the quantity /„ = p u /(p 7 + p v ). The 
effect of free streaming, even if more difficult to express 
analytically, is also mainly dependent on the value of f v 
[7(Tj . If we consider the case where the three standard 
neutrinos are the only contribution to the radiation en- 
ergy density other than photons, but we allow for the 
possibility of a non-vanishing chemical potential or for a 
different T°/T° ratio, we see again that the changes in 
the shape of the CMB anisotropy spectrum depend only 
on g eS as whole, as long as eq. @ retains its validity, i.e., 
as long as neutrinos are in their ultra-relativistic regime. 
Even considering the presence of some additional relic 
particle x does not seem to change this picture. Sup- 
posing that the other ultrarelativistic particles behave as 
neutrinos for what concerns the effects under considera- 
tion, we can argue that in all cases the relevant quantity 
is f x = p x j '/9 ra d, so that we are again lead to the con- 
clusions that g eS is the only relevant parameter. This 
means, for example, that in the case of neutrinos with 
mass less than ~ 0.3 eV, so that they stay ultrarelativis- 
tic until the time of last scattering and for some time af- 
ter, the effect on CMB perturbations is exactly the same 
of massless neutrinos, and every change in their temper- 
ature or chemical potential, as even the presence of an 
additional, sterile neutrino, is absorbed in g oS (moreover, 
we don't have obviously any possibility to extract infor- 
mation about their mass). 

However the picture changes when considering neutri- 
nos (or other relic particles) that go out from the ultra- 
relativistic regime before matter-radiation decoupling. If 
the neutrino mass is larger than ~ 0.3 eV, the effect of 
its finite mass is felt by the perturbations that enter the 
horizon after neutrinos have gone out from the ultrarela- 
tivistic regime, because from some point on the evolution 
the energy density will be given no more by the approx- 
imate formula 10, but instead by eq. 10) that contains 
a dependence on mass through the term (3. A side effect 
of this is that it will not be possible to single out the de- 
pendence of p from T and £ as an overall factor, so that 
these two contributions become distinguishable. 

Let us make this more clear with an example. Consider 
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a gravitational wave entering the horizon after neutrinos 
became non-relativistic, but before matter-radiation de- 
coupling. This wave will be absorbed, according to |69j |. 
proportionally to p v . On the other hand the free stream- 
ing length of neutrinos will vary according to the velocity 
dispersion < v 2 > |UE3|- The key point is that, for a 
gas of non-relativistic particles, p v and < v 2 > will de- 
pend on T v and in different ways, so that measuring 
independently the absorption factor and the free stream- 
ing length, it would be possible at least in principle to 
obtain the values of T v and without any ambiguity 
left. 

What we have just said is even more true with respect 
to the LSS data, since even neutrinos with mass greater 
than 10~ 3 eV are in their non-relativistic regime during 
the late stages of the process of structure formation. We 
conclude then by stressing that one should be careful 
when parameterizing the lepton asymmetry by means of 
an effective number of degrees of freedom. 



V. METHOD 

We used the CMBFAST code , modified as described 
in Ref. |46j | in order to account for a non vanishing 
chemical potential of neutrinos, to compute the temper- 
ature (TT) and polarization (TE) CMB spectra for dif- 
ferent combinations of the cosmological parameters. As 
a first step, we added three more parameters, namely 
the effective number of additional relativistic species 
A7V° t * crs , the neutrino degeneracy parameter £ (both de- 
fined in Sec. and the neutrino physical energy density 
uj v = n„h 2 to the standard six-parameters ACDM model 
that accounts in a remarkably good way for the WMAP 
data. As anticipated above, AJV^ crs accounts only for 
the extra energy density due to the presence of additional 
relic relativistic particles other than the three Standard 
Model neutrinos. We shall refer to the (cj„, £, AN^ hcrs ) 
subspace as the "neutrino sector" of the parameter space 
(although, as we have just noticed, AA^° t * crs does not 
refer directly to neutrinos). 

With the above mentioned choice of the parameters we 
can make a consistency check to our results, by verifying 
that imposing the priors £ = 0, u v — and AiV^ ers = 0, 
we obtain results that are in agreement with the ones of 
the WMAP collaboration. Moreover, by choosing a suffi- 
ciently wide range for the variation of the three additional 
parameters, we can check how much their introduction 
affects the estimation of the best-fit values of the core 
parameters. Thus we choose to use the following param- 
eters: the physical baryon density u>b = flbh 2 , the total 
density of non relativistic matter u) m = (fib + ^cdm)^ 2 , 
the scalar spectral index n, the optical depth to reioniza- 
tion t, the overall normalization of the CMB spectrum 
A, the physical neutrino density u) v = fl^h 2 , the neutrino 
degeneracy parameter £ and the extra energy density in 
non-standard relics AA/^L, rs . We will be considering the 
scenario in which the three standard model neutrinos 



have all the same mass and chemical potential. We take 
the chemical potential to be positive (this corresponds to 
an excess of neutrinos over antineutrinos), but since the 
effects on the CMB do not depend on the sign of £, we 
quote the limits that we obtain in terms of its absolute 
value. We do not include as a free parameter the Hubble 
constant Hq , whose degeneracy with the effective number 
of relativistic degrees of freedom and with the neutrino 
mass has been studied in previous works [f|. Instead we 
decided, according to the recent measurements of Hubble 
Space Telescope (HST) Key Project |x3 , to assume that 
h = 0.72. Moreover, we restrict ourselves to the case of 
a flat Universe, so that the density parameter of the cos- 
mological constant fl\ is equal to 1 — (ui m +u v )/h 2 . We 
are thus dealing with a 8-dimensional parameter space. 

Let us discuss in a bit more detail the way we deal with 
priors in the neutrino sector of parameter space, i.e., with 
information coming from other observations, and in par- 
ticular from BBN. As we have stressed in section ITT1 the 
standard BBN scenario, together with the equalization of 
chemical potentials, constraints the neutrino degeneracy 
parameter to values lower than the ones considered in this 
paper; on the other hand this conclusion possibly does 
not hold in non-standard scenarios where additional rel- 
ativistic relics are present. However, even non-standard 
scenarios of this kind usually single out some preferred 
region in parameter space. At the present, several non- 
standard scenarios that can account for the observed He- 
lium abundance exist (see for example Refs. and |l2() 
so that we adopt a conservative approach, and choose not 
to impose any prior on the neutrino sector, other than 
the ones that emerge "naturally" as a consequence of our 
choice of parameters. Anyway, this does not preclude the 
possibility of successively using the BBN information: in 
fact once the likelihood function in the neutrino sector 
has been calculated, it can be convolved with the rele- 
vant priors coming from non-standard BBN scenarios. 

We span the following region in parameter space: 
0.020 < uj b < 0.028, 0.10 < u m < 0.18, 0.9 < n < 1.10, 
< r < 0.3, 0.70 < A < 1.10, < u> v < 0.30, 
0<|CI<2.0,0< AiV o f hcrs < 2.0. We shall call this our 
"(5+3) parameter space". In order to obtain the likeli- 
hood function C(ui bl u> m , n, r, A, uj u , £, AN°f hcl . s ) in this 
region, we sample it over a grid consisting of 5 equally 
spaced points in each dimension. For each point on our 
grid, corresponding to a combination of the parameters, 
we compute the likelihood relative to the TT [l9j and TE 
pfjj angular power spectrum observed by WMAP, u sing 
the software developed by the WMAP collaboration |2jj 
and kindly made publicly available at their website 1 . To 
obtain the likelihood function for a single parameter, we 
should marginalize over the remaining ones. However for 
simplicity we approximate the multi-dimensional integra- 
tion required for the marginalization procedure with a 
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maximization of the likelihood, as it is a common usage 
in this kind of likelihood analysis. This approximation 
relies on the fact that the likelihood for cosmological pa- 
rameters is supposed to have a gaussian shape (at least 
in the vicinity of its maximum) and that integration and 
maximization are known to be equivalent for a multivari- 
ate Gaussian distribution. 



{HI together with the following relations El E3 : 



,2 _ m " Fi £v) 

Uul% ~^ 93.5 eV ' 

V 

_ n v — rip 



(20) 



According to Bayes' theorem, in order to interpret the 
likelihood functions as probability densities, they to be 
inverted through a convolution with the relevant priors, 
representing our knowledge and assumptions on the pa- 
rameters we want to constrain. Here we shall assume 
uniform priors, i.e. we will assume that all values of the 
parameters are equally probable. 



For each of the core parameters, we quote the max- 
imum likelihood value (which we shall refer to also as 
the "best-fit" value) over the grid and the expectation 
value over the marginalized distribution function. We 
quote also the best chi square value Xo ( we recall that 
X 2 = —2 ln£) divided by the number of degrees freedom, 
that is equal to the number of data (for WMAP, this is 
1348) minus the number of parameters. For what con- 
cerns the parameters of the neutrino sector, we quote the 
maximum likelihood and the expectation value as well, 
and in addition we report a 2a confidence interval. Us- 
ing a Bayesian approach, we define the 95% confidence 
limits as the values at which the marginalized likelihood 
is equal to exp[— (xo ~ 4)/2], i.e., the values at which 
the likelihood is reduced by a factor exp(2) with respect 
to its maximum value 2 . There is one exception to this 
procedure, namely, when the maximum likelihood value 
for a parameter that is positively defined (such as oj v 
or the absolute value of £), let us call it 6, is equal to 
zero. In this case, instead than computing the expec- 
tation value, we just give an upper bound. In order to 
do this, we compute the cumulative distribution function 

= ( So £(9)d0) I (ja° C(9)d0) and quote as the up- 
per limit at the 95% confidence level the value of 9 at 
which C{9) = 0.95. 



12C(3) 



E(^)(f) 



(21) 



where 



no 



3C(3) 



k=l 



3C(3) 



1 ^ 2 °° 



k=l 



fc 3 



(22) 



VI. RESULTS AND DISCUSSION 

We start our analysis by looking at the effect of the 
introduction of the additional parameters to the estima- 
tion of the core parameters (cob, ui m , n, r, A). First of 
all, we check that imposing the priors £ = 0, lo u — and 
^-^othcrs = our results are in good agreement with the 
ones of the WMAP collaboration (we should refer to the 
values quoted in Table I of Ref . pj ) . The mean and max- 
imum likelihood values that we obtain for each parameter 
are summarized in the second and third column of Table 
|U We see that in all cases our results lie within the 68% 
confidence interval of WMAP expected values. Then we 
remove the prior on u> v , while still retaining the ones on 
£ and AiV^jj . The maximum likelihood model has still 
u) v = 0. The best- fit values of the core parameters are 
left unchanged, and the same happens for the best-fit x 2 , 
thus suggesting that a non-zero lo v is not required in or- 
der to improve the goodness of fit. The results for the 
core parameters are summarized in the fourth and fifth 
column of Table [I] 



Once we have obtained constraints on uj u , £ and 
ANqR , we translate them to limits on the neutrino 
mass m v , the lepton asymmetry L and the extra number 
of effective relativistic species AiV eff , using eqs. l(T%|l and 



2 The 95% confidence level defined in this way is not in general 
equal to the 2a region, defined computing the variance of the 
probability distribution. However, the two are equal for a gaus- 
sian probability density. As we shall see, almost all the marginal- 
ized distribution have a nearly gaussian shape. When it is not 
so, we shall point this out. 



Finally, we compute the likelihood over our whole 
parameter space. The results for the core param- 
eters are summarized in the sixth and seventh col- 
umn of Table [I] The maximum likelihood model 
over the grid has (ui b , uj m , n, t, A, uj u , £,,AN^ hcls ) = 
(0.022, 0.14, 0.95, 0.075, 0.7, 0, 0.5, 0). We see that this 
time, the best fit values for the five core parameters 
are slightly changed with respect to the standard case. 
The changes in uj m and n could seem strange at a first 
sight, since intuitively one would expect the opposite be- 
haviour, i.e., a change to larger values for both, because a 
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TABLE I: Core Parameters 



Priors: f = 0, Q. u = 0, AJV^f, erB = Priors: £ = 0, AJV°^ era = No Neutrino Priors 

Maximum Maximum Maximum 

Parameter Mean Likelihood Mean Likelihood Mean Likelihood 



Baryon density,^ 


0.024 


0.024 


0.024 


0.024 


0.023 


0.022 


Matter density, u; m 


0.15 


0.16 


0.15 


0.16 


0.14 


0.14 


Hubble constant a ,ft, 


0.72 


0.72 


0.72 


0.72 


0.72 


0.72 


Spectral index, n 


1.00 


1.00 


1.00 


1.00 


0.98 


0.95 


Optical depth, r 


0.13 


0.075 


0.13 


0.075 


0.12 


0.075 


Amplitude, A 


0.8 


0.8 


0.8 


0.8 


0.8 


0.7 


x 2 /" 




1437/1343 




1437/1342 




1431/1340 



"The value of the Hubble constant is kept fixed to h = 0.72. 



6000 




1 1 1 ' ■ — ■ ■ — 1 

1 10 100 1000 

Multipole Moment I 



FIG. 1: Comparison between the best-fit power spectrum ob- 
tained assuming the priors £ = 0, AA^ t ff hcrs = (solid line) and 
without such prior (dashed line). The points are the WMAP 
data on the temperature angular power spectrum [Tg|. 



larger oj m could keep the time of matter-radiation equal- 
ity, while a larger n would increase the power on small 
scales thus leaving more room for neutrino free stream- 
ing. This is because the goodness of fit of a particular 
model with respect of the WMAP data is mainly deter- 
mined by its ability to fit the first and second peak. In- 
creasing together ui m , n, £ and AiV° t ^ ers would increase 
the height of the first peak that can be then lowered 
back by decreasing the overall amplitude A. We show in 
fig. ^a comparison between the best-fit spectrum in the 
(£ = 0, AJV^ lcrs = 0) subspace with the best-fit spec- 
trum on the whole space. 

Now let us turn our attention to the neutrino sec- 
tor of parameter space. The best-fit model over the 
(£ = 0, AA^ lers = 0) subspace of the grid has w v = 0, 
and the x 2 changes from 1437 to 1541 when going from 
u) v = to the next value in our grid, us u — 0.075. We 
can compute an upper bound for uj v , but since the region 
in which C(lo v ) significantly differs from zero is all com- 
prised between the first two values in our grid {0, 0.075}, 
the result is rather dependent from the particular inter- 
polation scheme we choose. Using a simple, first order 
interpolation scheme, we find the bound uj v < 0.0045 



(95% CL), corresponding to Tfiu <C 0.14 eV, while using 
higher order interpolation schemes the bound weakens up 
to w v < 0.015 (m„ < 0.47 eV). This result should then 
be taken with caution and we shall simply consider it as 
an indication that, although we are using a grid-based 
method with a rather wide grid spacing instead than the 
more sophisticated Markov Chain Monte Carlo (MCMC) 
method [73l F3 . Fij , we basically obtain the same results 
of the WMAP collaboration, namely, w„ < 0.0072 Q, 
when imposing the priors £ = 0, AiV° t * crs = 0. 

We make a second check by imposing that w v — 0, £ = 
and computing the 95% confidence region for AiV^ crs . 
We find that < AiV^ crs < 1.4. Since the degeneracy 
parameter is vanishing, the same limit applies to AA?" eff = 
A-^otnors- This is quite in agreement with the results 
quoted in Ref. although it is more restrictive. This is 
probably due to the fact that we are imposing a stronger 
prior on h, keeping it constant and equal to 0.72. This is 
confirmed by a visual inspection of fig. 2 of Ref. . 

The best-fit value for neutrino density over the whole 
parameter space is still w„ = 0, but this time % 2 changes 
from 1431 to 1441 as ui v goes from to 0.075, so that 
the probability density spreads out to higher values of 
uj v with respect to the preceding case. The result is 
that the upper bound raises up to lo v < 0.044, quite in- 
dependently from the interpolation scheme used. This 
is probably related to the already observed trend for 
which, when the energy density of relativistic relics is 
increased, the pos sibility for larger neutrino masses re- 
opens IliEIIllzl 

The maximum likelihood value for the degeneracy pa- 
rameter is |£| = 0.5, while the expectation value over 
the distribution function is |£| = 0.56 (corresponding to 
\L\ = 0.43). At the 2a level, the degeneracy parameter is 
constrained in the range < |£| < 1.07. This corresponds 
to < \L\ < 0.9. For what concerns the additional num- 
ber of relativistic relics, the maximum likelihood model 
has AJV"^^ = 0, and the expectation value over the 
marginalized probability function is AA r ° t f J lcrs = 0.3. The 
95% confidence region is -0.7 < AiV ™ crs < 1.3. This 
opening towards smaller, negative values of AiV° t f ^ lors can 
be ascribed to the fact that such values produce a lower- 
ing of the acoustic peak, that can be compensated by a 
larger degeneracy parameter. The quoted bounds on u v 
and £ translate to the following bound on the neutrino 
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0.2 0.4 0.6 

m v (eV) 




FIG. 2: Likelihood functions for ui v , |£|, AN^f hcls and for the derived parameters \L\, ra„ and AN cB , obtained as a result of 
the analysis in (5+3) parameter space. The dotted lines bound the 95% confidence interval. The functions are normalized so 
that their integral is equal to unity. 



mass: m„ < 1.2 eV (95% CL). In fig. [2] we show the 
likelihood functions, while in Table ITT1 we summarize our 
results for the basic and derived parameters describing 
the neutrino sector. 

We remark that, although the maximum likelihood 
model over the whole grid has uj v = 0, this is not in con- 
tradiction with our choice of considering £ and AiV^ s 
as independent parameters, in spite of the fact that in 
this limit they should be degenerate. The basic reason is 
that, as can be seen from the likelihood curves, models 
with io v > can be statistically significant. For these 
models, AN° S and AN^ heTB are not exactly degenerate. 

In order to better study the partial degeneracy be- 
tween |£| and AN^ hers , and then to understand how the 
value of AJV^ crs affects the estimation of the degener- 
acy parameter, we compute the likelihood curve for the 
degeneracy parameter for particular values of AN^ hcrs . 
The results are shown in Table ITTTl . From this table, a 
quite evident trend appears, namely that for large values 
of AN^ hcrs , smaller values of |£| are preferred, and vicev- 
ersa. As already noticed, this is probably related to the 
fact that when AiV^ff is increased, it remains less room 
for the extra energy density of neutrinos coming from the 
non-vanishing degeneracy parameter. It is worth noting 
that, for AN^ hcls ~ 0, the case £ = lies outside the 
95% confidence region. We stress the fact that, accord- 
ing to theoretical predictions, in models of degenerate 
BBN with "3+1" neutrino mixing, if chemical potentials 
are large (£ > 0.05), the production of sterile neutrinos is 
suppressed, effectively resulting in AN^ OIS = 0, 0] . 



TABLE II: Neutrino Sector 





95% Confidence 


Parameter 


Interval 


Physical neutrino density, uj v 


< 0.044 


Degeneracy parameter, |£[ 


60 +0 - 50 


Neutrino mass in eV, m„ 


< 1.2 


Lcpton asymmetry, \L\ 


46+ 43 

u -^° — 0.46 


Effective number of additional 


0.30 ± 1.0 


rclativistic relics, A]V° t f [ ors 




Effective number of additional 




rclativistic relics, AiV cft 




TABLE III: Correlation between £ and AiV*^ ers 


A /V cff 
"-"others 


i 

(95% Confidence Interval) 





0.65 ± 0.58 


0.5 


42+ 58 


1.0 


o-is±S:?| 


1.5 


< 0.53 


2.0 


< 0.29 



VII. CONCLUSIONS AND PERSPECTIVES 

In this paper, we have studied the possibility to con- 
straint the lepton asymmetry of the Universe, the sum 
of neutrino masses, and the energy density of relativis- 
tic particles using the WMAP data, in the framework of 
an extended flat ACDM model. Despite the fact that 
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the current amount of cosmological data can be rather 
coherently explained by the standard picture with three 
thermally distributed neutrinos, vanishing lepton asym- 
metry and no additional particle species, nevertheless we 
think that it is useful to explore how non-standard sce- 
narios are constrained by the cosmological observables. 
We have concentrated our attention to models with a 
(eventually large) net lepton asymmetry (corresponding 
to a non-zero degeneracy parameter for neutrinos). Such 
models are motivated in the framework of extensions to 
the standard model of particle physics, and can possibly 
explain the observed amount of baryon asymmetry in the 
Universe. Having in mind this, we have also included the 
energy density of relativistic species as an independent 
parameter. In this last aspect, our approach differs from 
previous ones, where the two parameters where consid- 
ered degenerate. We have remarked that, although an 
approximate degeneracy between the two exists, it could 
be broken by finite mass effect, especially in the case of 
neutrino masses saturating the tritium beta decay bound. 

When considering perfect lepton symmetry , our re- 
sults are in agreement with previous ones. In the more 
general case, we have found that, at the 2a level the 
bounds on the degeneracy parameter and lepton asym- 
metry are respectively < |£| < 1.1 and < \L\ < 0.9. 
The effective number of additional relativistic species (ex- 
cluding the contribution from the non standard thermal 
distribution of neutrinos) is bounded as follows (95% 
CL): -0.7 < AA^ lers < 1.3. Including also neu- 
trinos, this reads -0.45 < AN cS < 2.10. This limit 
is much more restrictive than the ones found in similar 
analysis 0,^3- This is probably due to the fact that we 
assume a very strong prior on the Hubble parameter, fix- 
ing h — 0.72. The physical explanation is that the later 
matter-radiation equality due to AN cS > can be com- 
pensated by making u> m = fl m h 2 larger, and viceversa. 
This gives rise to a partial degeneracy between AA off 
and h, thus making the costraints on both parameters 
looser unless some external prior is imposed to break the 
degeneracy. 

We also find that the data are compatible with u v and 
vn v equal to 0, with upper bounds (95% CL) u) v < 0.044 
and m v < 1.2 eV. This bounds are larger than the ones 
usually found, and this is probably due on one hand to 
the presence of a larger energy density of UR particles, 
and on the other hand to the wide grid spacing we have 
used. 



The usual scenario, with \L\ = and AN cB — 0, is 
then compatible with WMAP data at the 2a level; how- 
ever the likelihood curves show that alternative scenarios 
with £ ~ 0.6 and AN cS ~ 0.7 have a larger likelihood 
with respect to the data. In effect, the standard sce- 
nario lies outside the la confidence region. Even if this 
is not enough to definitely claim evidence, in the CMB 
anisotropy spectrum, of exotic physics, we think that it 
is however interesting that non-standard models are not 
ruled out but actually preferred by the WMAP data. 

We have also studied how the results on the lepton 
asymmetry can change when more precise information on 
the energy density of relativistic particles is given. We 
have shown that, the smaller is the extra energy density, 
the larger is the allowed lepton asymmetry. In particu- 
lar, for models with vanishing AN^ beis , perfect lepton 
symmetry is ruled out at the 2a level. This is probably 
due to the approximate degeneracy between AN eS and 
£. The issue of the exact extent of this degeneracy is still 
open, and we think that it deserves a deeper attention. 
It would be desirable to investigate if future precision 
CMB experiments, and in particular the PLANCK mis- 
sion, can clearly disentangle the two parameters. Even 
more promising to this purpose is the power spectrum of 
LSS, as we have stressed in Sec. IIVI 

The results presented in this paper have been derived 
assuming that the three neutrino families have all the 
same mass and chemical potential, owing to the struc- 
ture of the CMBFast code used to generate the theoretical 
power spectra. It is yet to investigate the role of a non- 
uniform distribution of the lepton asymmetry between 
different families. In a similar way we did not investigate 
models in which chemical equilibrium between neutrinos 
and antincutrinos does not hold, implying £„ =/= —£5. 

We conclude then that WMAP data still cannot ex- 
clude the presence of non-standard physics in the early 
evolution of the Universe. In particular, they do not 
exclude the presence of a large neutrino asymmetry, and 
consequently they do not rule out exotic leptogenesis sce- 
nario where a large lepton number is produced. 
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